Convergence method for calculating solutions to the 3D invariant 
embedding integro-differential equations describing electron 

transport processes. 

C. Figueroa, H. Brizuela, and S. P. Heluani 

Laboratorio de Fisica del Sdlido, Dpto. de Fisica, 
Facultad de Ciencias Exactas y Tecnologia, 
Universidad Nacional de Tucumdn, Argentina. 
(Dated: June 22, 2010) 

Abstract 

The electron and photon transport processes in spectroscopy techniques described by the invari- 
ant embedding theory is here revisited. We report a convergence method to obtain closed analytical 
solutions to the 3D integro-differential equations. This method was successfully used in calculating 
the dependence of the electron backscattered fraction on the atomic number and on the energy. 
Also the fraction of absorbed electron as a function of incident angles was calculated. Using a states 
ladder model for the electron energies, this method provides a tool for testing physical parameters 
involved in the transport theory, such as the elastic and inelastic cross sections. The outstanding 
feature of the invariant embedding differential equations of considering observable quantities (such 
as the emergent flux of particles) as independent variables makes them a suitable tool to describe 
experimental situations. 
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I. INTRODUCTION 



Advances in characterization techniques apphcable to the complex structures of materials 
such as nanorods, thin films, nanofiltration membranes etc., will advance the understand- 
ing of the relationship between structure and properties in the search for technological 
applications. In the field of structural and chemical characterization methods, the elec- 
tron and photon transport has several applications including Electron Probe Microanaly- 
sis (EPMA), X-ray Fluorescence Analysis, Electron-Beam- Induced-Current, Auger Electron 
Spectroscopy, Electron Energy Loss Spectroscopy, X-ray Photoelectron Spectroscopy (XPS), 
and Secondary Electron Emission. All these characterization techniques are based on the 
strong interaction of electrons and/or photons with matter and in the emissions or 'signals' 
that they produce. The different applications are defined according to which of these signals 
are detected and how they are analyzed. Particularly, in quantitative microanalysis, the 
energies and intensities of the characteristic x-ray lines from an irradiated sample are used 
to determine the elemental composition or dimension of small samples and also the films 
thicknesses. To quantify it the 'signal' must be corrected using some transport model and 
also some 'correction model'. The loss of ionization due to backscattered electrons is one 
of the aspects that the correction models have to take into account and an example of the 
complexity of the problem. Therefore, most of the approaches which form the basis of the 
correction procedures commercially available in spectroscopy techniques, are still empirical 
or semi-empirical ones^. 

Theoretical and experimental studies of the backscattered electrons are valuable in sev- 
eral techniques, such as reflection Electron Energy Loss Spectroscopy, Particle Induced X- 
ray Emission and Elastic Peak Electron Spectroscopy. Also their applications in Scanning 
Electron Microscopy (SEM) and in film thickness determinations are of great importance. 
In particular, SEM represents a high-performance method for investigating structures and 
devices in the nanometer scale^"-. 

To deduce the presence of a particular chemical element, or to measure the thickness of a 
film or the size of a particle from spectroscopy measurements, it is necessary to describe the 
emergent flux of particles (an observable quantity) in terms of the interactions produced by 
the injected flux during its trajectory inside a solid sample. Several mathematical models 
were developed to describe the evolution of the electrons or photons in a solid medium. 



2 



The application of the Boltzmann equation, or similar transport equations, presents serious 
difficulties in solving the resulting integro-differential equations and little progress has been 
made in the solution of 2D and 3D transport problems. In addition more difficulties arise 
in the practical applications of theoretical results to different techniques^"®. For example, 
Bernasconi et al." consider the hot-electron transport in the conduction band of a thin in- 
sulating film in terms of the energy-dependent elastic and inelastic scattering rates. The 
authors obtain an exact formal solution for a system of integro-differential Boltzmann equa- 
tions only for the zero-energy-loss transport problem. Werner- considered the backscattering 
of medium energy electrons from solid surfaces analyzing a linearized Boltzmann-like kinetic 
equation. Then, an approximate solution is implemented in a Monte Carlo scheme. 

The Boltzmann transport equation and Monte Carlo calculations (which also make use 
of the Boltzmann transport equation to simulate the track of photons or electrons) are still 
the most used methods for the interpretation of the signals in spectroscopic and characteri- 
zation techniques. However, they are uneconomical methods, since much of the information 
obtained about the internal fluxes is quite useless to the experimenter. More rigorous treat- 
ments refiexing the complexities of the problems are the semiclassical treatments and Fresnel 
optics equations^"~ and quantum treatments'^. 

In recent articles the so-called Invariant embedding Approach to Microanalysis (HAM) 
was developed with the aim of obtaining analytical solutions to the transport prob- 
lems associated to microanalysis techniques. The HAM is free of the drawbacks of the 
empirical methods and the mathematical difficulties to solve the Boltzmann equations 
(see ref.-). Using a simple ID model, expressions for the detected X-rays characteris- 
tic intensities had shown a good performance in the interpretation of experimental data.— 

The invariant embedding method is a mathematical technique based directly upon the 
physical processes. In this frame, the traditional linear functional equations, subject to 
boundary conditions, are transformed into nonlinear functional equations subject only to 
initial conditions in space and time coordinates. Although analytical approaches to solve 
invariant embedding differential equations were reported in different physics contexts, the 
method has, as yet, had relatively little impact in spectroscopy techniques^"—. 
The purpose of this article is to report a convergence method for deriving an analytical 
solution to the 3D invariant embedding equations associated to a states ladder model for 
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electrons transport processes within a solid medium. Considering the convergence of an ap- 
proximate solution to the exact solution in a medium of infinite size, the distributions of the 
backscattered and the absorbed electron fractions are obtained. In section II a brief descrip- 
tion of the geometrical configuration and the parameters used in the theory is presented. 
Then, analytical expressions for the electron fluxes are obtained as functions of physical and 
experimental parameters such as the ratio between the elastic and inelastic cross sections, 
the incident energy and the incident and emergent angles. Finally, the theoretical results 
are compared with experimental data. 

II. METHOD OF CONVERGENCE. 
A. Parameters of the theory. 

Considering a typical experimental design in EPMA, one of the advantages of the in- 
variant embedding method is to take the sample thickness as the integration variable in the 
differential transport equations (instead of the coordinate measured from the sample surface, 
as usual in the traditional Boltzmann equation treatment). Then, it is helpful to keep in 
mind the fact that the theoretical expressions for the different fluxes of particles (such as 
backscattered electrons) are written as a function of the coordinate measured from the right 
end of the sample which is also the thickness of the sample, see flgure lb. 

Let us consider an usual experimental conflguration in quantitative electron spectroscopy 
methods, as sketched in Figure 1-b, where a solid sample of a flnite thickness r is irradiated 
from the left by a beam of electrons with energy Eq (this geometrical conflguration could 
be applied to other similar techniques and the electrons could be photons). Some of the 
impinging electrons of the primary beam could be backscattered, some others could become 
deactivated (i.e. they become unable to produce ionizations or characteristic x-rays inside 
the sample) and some others could escape from the other side (in the case of fllm samples or 
light particles). In passing through the material the electrons undergo elastic and inelastic 
collisions. Both interactions are described by cross sections that rigorously are energy and 
angle dependent. In this work we consider the usual axial symmetry and only azimuthal 
integrated cross sections is considered. 

To clarify the notation, the diagrams in Figure la shows the angles {aj and /3j), which 
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are the incident and emerging polar angles for a single interaction, and the index j = 0, 1, .. 
labels the distinct energy states values. The polar angles are measured with respect to an 
axis perpendicular to the sample surface. 

For isotropic elastic cross sections, the integral with respect to the azimuthal angle leads 

to, 

dap- = 2tt k sinfSjdPj. 

Where k is a constant. The integrated cross section is: cxj = J dap. = Air k. 

To simplify the notation we shall call cos/3j = (3j. Considering only the energy level Eq, 

dcTpo = ao/2 d(3o. (1) 

Let us define the ratio of the elastic and the inelastic cross sections as X{E) = a{E) / s{E) 
, where a{E) is an estimation of the elastic cross section per path length unit and s{E) is the 
probability per length unit that the electron be inelastically scattered. The isotropy of the 
elastic cross section makes the expressions angle-independent. The ratio A has the advantage 
that it could be considered as energy independent, making the invariant embedding equations 
more simple to manage. 

In this work, the states ladder model described in refi^ will be used considering five 
energy states. The energy range of active electrons (i.e. capable of produce ionization) is 
divided into intervals or steps in a ladder. The transition (in the ladder) from one step to 
the following is described as follows: after an inelastic collision the electrons change their 
energies by a discrete amount 5E^ always degrading its original value, which implies to 
consider the problem of multiple inelastic collisions as reduced to a problem of only a single 
effective collision. 

B. Invariant embedding equations 

As usual in lEAM if it is possible to add a new layer of infinitesimal thickness dr on 
the surface of the sample, then we may consider the interval (0, r) as a sub-sample lying 
in the new (0, r + dr) sample^ii^ii^i^S. Let us introduce the matrix R = {Ra^^p,.) where 
Raj,i3k is the probability that an electron that impinges on the sample with polar angle 
aj and with energy Ej leaves the sample with an angle (3 and energy Ek- Considering 
only one energy state, figure l.b shows all possible trajectories which make contribution 
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to the electron backscattered fraction from a sample of thickness r + dr. Consider as an 
illustration the probabilities of the occurrence of the paths associated to the first and second 
diagrams in figure 1. The probability of occurrence of the first path is the probability 
that an impinging electron having an energy Eq and cosine of the polar incident angle 
ao passes across the interval (r + dr, r) without suffering elastic collision neither inelastic 
collision, it is (1 — sodT/ao){l — (Todr/ao) multiplied by the probability that the electron 
will be backscattered from a "sub-sample" of size r without loss its energy, it is -Rao,/3o 
by the probability that leave the sample at t + dr with an emerging cosine Po without 
suffering collision in the interval (r + dr^r). The probability of the path is written as 
(1 — soc?r/ao)^(l crodT/(^oYRao,i3o- Using the same procedure, the probabilities of the 
other four paths are written as: 

Path 2 (1 - so^Jda^,,So^^Rso,Po,a ' So%){l - 

3 (1 - s,^){l - ao^jRao,eJl - s,%)da,,,p,% 

4 (1 - .o^)(l - ao^)i?«„,,Jl - 

5 (1 - 5o|^)^^^^ao,/3oS• 
It is important to emphasize the fact that the details of the particle behavior in (0, r) are 

of no interest. As usual, to obtain i?ao,/3o('^ + dr) all probabilities of the paths illustrated 
in Fig. 1 must be summed taken into account that the probabilities of the trajectories 2, 
3 and 4 must be integrated over the cosines eo and ip^. Then letting dr — )■ to yield to 
a differential equation for i?ao,,9o(''")^^- ^'^Y other possible paths have a contributions of the 
order of {drY and their contribution will vanish in this procedure. The differential equation 
for Rag^pg, in its simplest form for one level model is written as. 



dr ao '^0 

do'eo,^oRao,eo {'t)R'Po,Pq ij) 



, daeo,l3oRaQ,e„ir) w Sq + (To Sq + CTq , , s 

+ / RaoM^)[ n J- (2) 

eo Po tto 



The integro differential Eq. 2 do not have an exact solutions. However, using a a 
convergence method it is possible to obtain a practical solution to be used in spectroscopy 
measurements. 
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Suppose T oo for the trajectories in Fig.l that involve at least one coUision in the 
differential layer dr. These are the trajectories 2, 3, 4 and 5 in Fig.l. Let us introduce 
the function r{aQ, Pq, oo) which denote the probability of occurrence of these trajectories so 
that, 



r(«o,^o,oo) = ^^+ / 
"0 Js, 



So Q^O 
dcTeo ,yo Rgp ,60 ( OO ) -Ryp , ft ('^ ) 

eo 



_^ f ^^eo,;3o-Rao,eo(oo) ^g-j 

If we now replace the first four terms in Eq.2 by the later expression, an approximate solution 
for RaQ,/3Q{T) can be written as. 

We shall now proceed to derive r(ao, /3o, c>o) which makes the approximate solution 
i?ao,/3o Eq.4, to converge to the exact solution of Eq.2 for large values of r. 



C. Convergence of the solution 

Let us now proceed to replace i?o,(,^^o(oo) from Eqs 4. in Eq. 3 and then we proceed to 
solve the following integrals, 

r ao,/3o,oo>i/3o = ^r— X {1+ / 77-— ttt-— — v + 

'-^ deor{ao,/3o, 00)010 



r f 

Jo 



j 

Jo 



(cKo + eo)(so + (To) 



I / deod(/?or(ao,^o,oo)^Q;o(/?o/3o 
' ' {ao + eo){<po + l3o){so + ao)'^' ^ ^ 

where the cosines Po, ^0 a-nd (po are represented in Fig. 1-a. In Eq. 5 we make use of the 

fact that the integration variables arc the director cosines of the emerging or of the incident 
angles, which by symmetry can be exchanged. 

Exact evaluation of the integrals in Eq. 5 is a difficult task, because it is necessary to know ex- 
plicitly the angular dependence of r{ao, Po, 00) x ao- To proceed, taken into account that the 
unknown function is a continuous function of ckq, Po for all angle, we consider that the value 
of r{ao, Po-i 00) X Uo in the integrands of Eq. 5 is nearly independent of the director cosines 



and equal to a constant which must be estimated using physical arguments. One method 
could be to use the mean-value theorem. However, in this work we consider an effective 
value r* that could be calculated from the following procedure: substitute r(ao, /3o, oo) x ao 
by the unknown constant r* in the arguments of the integrals in Eq.5 and then evaluate the 
right hand member. The result for r(ao) Po-, oo) is, 

ao^(ao, /3o, oo) = (To/2 + ^ ^0*^0 . Ln{l + 

2(so + cToj Po 

2(so + (To) ao 2(so + o-q) 

Ln(l + -)Ln(l + l). (6) 

ao Po 

Now we make use of this later results to find the value of r* that satisfies the following 
conditions (considering the integrals in Eq. 5) , 



1 rl 



r(ao, /3o, oo)ao/3o'^aoc?/3o r* (3odaod(3Q 



lo Jo (ao + /3o)(so + cro) Jq Jq {ao + (3o){so + ao) 
Then, using (6) and (7) the final expression for -Rao,/3o(r = 00) in Eq.4 is, 

''-''•<~> = (A + l)(a„ + A) " f— 

' ' + l//3o)Ml + l/«o) 



(7) 



A 

+/3o'i:ri(l + l//3o) + Tiaol3oLn{l + l/«o)}. (8) 



where Fi = A + 2 — 2-\/A + 1, and C is a normalization constant which considers the conser- 
vation of the particle fluxes. 

Another useful function in characterization techniques is the fraction of absorbed elec- 
trons. Defining the matrix A = (Aaolr)), where Aa^ir) is the probability that an electron 
that impinges in the sample with cosine ao and with energy Eq becomes deactivated (i.e., it 
loses its capacity to produce ionization) inside the sample. This probability can be obtained 
using the convergence method outlined above. To do this, we have to construct diagrams 
similar to that shown in Fig. 1-b. (see RefiS for details). The resulting equation is: 

dAaoir) _ daao,l3oAl3o{T) , £0_ _^ /" SoRao,l3o{T) 



dr ao ao Jso 

eo Po ao 
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The solutions are: 

^ , , a(oo)a;o 

^ao(oo) = ^ ' ■ (10) 

So + 0^0 

An expression for a(oo) is easy to obtain using the same procedure as for r(oo) 

A., = X [1 + + ^Ln(l + l/ao)(l + T,)]. (11) 

D. Numerical evaluation 

On the basis of the procedure described above a set of closed expression for Raf^.p^ij), k = 
1, 2, 3, 4 and 5, is obtained using the states ladder model described in Ref. 14. The algebra 
is cumbersome but using a standard PC the calculation can be performed in a few minutes 
using any mathematical software. The method was applied to a number of elements of 
representative atomic numbers to illustrate its qualitative behavior. 

The numerical results depend on the approximations used to estimate the cross sections 
for the scattering processes involved. Although any approximation could be used in HAM 
equations, in this work the Rutherford elastic cross section is usect^^. The inelastic cross 
section is estimated from Bethe theory^^. Then, the ratio A could be expressed as: 

A = 3.96 10-20 i — \11A^L^_ 12 

where v is the minimum value of energy loss in each inelastic collision, / is the ionization 
energy of the atoms, p, N a and A are the density, the Avogadro constant and the atomic 
number respectively. 

The backscattered fraction 77, calculated as a function of the atomic number Z and the 
energy distribution for various materials are show in Figure 2a and 2b, respectively. Exper- 
imental data in figure 2a were taken from Database of electron- Solid Interactior?^ . Figure 
2a shows theoretical calculation considering i/ as a fitting parameter in the ratio A (or in 
the cross section expression). The best fit is obtained with v close to 4.25 xlO^^^KeV. 
The theoretical results lead to reasonably accurate agreement with the experimental data 
for the atomic number dependence. However, in the case of the energy spectrum for the 
backscattered electrons, our theoretical model predicts a maximum at the incident energy 
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EO. Actually, in experimental results the maximum in the spectrum occurs at energies 
slightly below EO (see Ref.— ). The discrepancy could be explained taking into account that 
our model considers the possibility of perfect elastic collisions, while actual collisions occur 
with certain amount of energy loses. 

Even though the limitations of the convergence approach to describe the maximum ex- 
pected in Fig. 3a, Eq. 2 was derived from first principles, straight from the model of the 
physical process. Therefore, it does not contain the strong physical hypothesis assumed in 
most of the models employed in different software packages for EPMA, namely, that the 
backscattered electron trajectories and their energy loss mechanisms inside the solid, are 
independent of the trajectories of the deactivated (absorbed) electrons. 

A qualitative behavior of the theoretical results for the fraction of absorbed electrons as 
a function of incident angle is presented in figure 3. In agreement with experimental results, 
the maximum occurs at normal incidence of the electron beam. 

III. SUMMARY 

The aim of this work is to report a new approach to solve the functional equations of 
IE that describes physical processes in spectroscopic techniques. With this approach it 
is possible to obtain, with relative facility, analytical expressions of accessible treatment 
which could be useful to the experimental investigators to interpret their results. These 
expressions also facilitate the evaluation of such parameters like cross sections, attenuation 
factors, etc. In the present work it we have dealt with the calculation of the absorbed and 
backscattered electronic fractions, improving previous models and thus allowing a closer 
approximation to the real phenomenon. This approach offers an different point of view for 
the study of the already mentioned physical parameters. At present there are calculations 
in progress of other parameters of interest in spectroscopy, such as the k reasons^ applying 
the convergence method. This study will allow an independent validation of the model and 
method previously described. 
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Figure captions 



Fig. 1: a) Angles defined for equation 3. b) Five trajectories wliicli make contributions to 
the backscattered electron fraction in a one-state model. Here a, ^ etc indicate the cosines 
of the corresponding angles. 

Fig. 2: a) Backscattered electrons fraction, 77, as a function of atomic number. Theoretical 
values, color points. Experimental values, crosses, b) Energy spectrum of backscattered 
electrons from samples of different atomic number Z. 

Fig. 3:The fraction of absorbed electrons as a function of incident angle. 
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FIG. 1: :a) Angles defined for equation 3. b) Five trajectories which make contributions to the 
backscattered electron fraction in a one-state model. Here a, (3 etc indicate the cosines of the 
corresponding angles. 
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FIG. 2: a) Backscattered electrons fraction, rj, as a function of atomic number. Theoretical 
values, color points. Experimental values, crosses. b)Energy spectrum of backscattered electrons 
from samples of different atomic number Z. 
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